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INTRODUCTION 


An  analytic  method  is  not  available  for  the  computation  of  stress 
intensity  factors  for  multiple-radial  cracks  in  a  thick-walled  cylinder.  The 
computation  must  depend  on  various  numerical  methods. 1“^  Due  to  increasingly 
successful  applications  of  the  finite  element  technique  in  structural 
analysis,  the  author  decided  to  use  higher  order  finite  elements  with  the  aid 
of  special  crack-tip  elements  to  study  multiply  cracked  cylinders.  Similar  to 
the  quarter-point  element  in  an  8-node  quadrilateral  element, ^>7  a  special 
crack-tip  element  was  developed^  for  a  12-node  quadrilateral,  isoparametric 
element.  Both  the  8-node  and  12-node  quadrilaterals  have  been  implemented  in 
the  popular  finite  element  computer  code  NASTRAN.**^  The  dummy  user  element 
facility  of  NASTRAN  is  used  for  the  implementation.  Another  finite  element 
computer  code  APES,l^  which  was  written  specifically  for  the  use  of  12-node 
quadrilateral,  isoparametric  elements,  has  also  been  used  for  fracture 
analysis  of  cracked  hollow  cylinders.  Quite  accurate  results  of  stress 
intensity  factors  have  been  obtained  using  either  NASTRAN  or  APES  for 
multiple-radial  cracks  emanating  from  the  bore  of  a  tube.!*  These  results  are 
in  good  agreement  with  results  reported  by  Tracy  12  using  the  method  of 
modified  mapping  collocation. 

To  increase  the  maximum  internal  pressure  a  cylinder  can  contain 
elastically  and  retard  the  growth  of  radial  cracks  near  the  bore,  it  is  a 
common  practice  to  introduce  compressive  residual  stress  near  the  bore  by  an 
autofrettage  process.  However,  the  residual  stress  in  the  cylinder  causes 

References  are  listed  at  the  end  of  this  report. 


1 


increased  difficulties  in  the  estimation  of  stress  intensity  factors.  One  of 
the  difficulties  is  that  results  of  residual  stress  distribution  vary  based  on 
different  assumptions  by  different  investigators.  Another  difficulty  is  the 
lack  of  an  initial  stress  analysis  capability  in  NASTRAN  and  APRS.  A  method 
is  developed  in  Reference  13  so  that  NASTRAN  or  APES  can  be  used  for  the 
computation  of  stress  intensity  factors  for  cracked  cylinders  with  residual 
stress  distribution  given  in  the  closed  form  expressions.^  The  finite 
element  results^  are  in  close  agreement  with  Parker's  results^  using 
modified  mapping  collocation. 

While  the  autofrettage  process  produces  favorable  compressive  residual 
stress  near  the  bore,  it  also  yields  a  tensile  residual  stress  near  the  outer 
cylindrical  surface.  The  sum  of  this  stress  and  the  tensile  stress  due  to  a 
bore  pressure  may  be  high  enough  to  cause  crack  initiation  and  propagation 
from  the  outer  surface  of  the  cylinder.  This  requires  the  computation  of 
stress  intensity  factors  for  externally  cracked  cylinders.  Several 
investigators  have  reported  their  results  on  this  subject. 12, 17-19 

In  order  to  reduce  repeated  finite  element  computations,  the  weight 
function  method^O  is  used  together  with  the  finite  element  method.  In  this 
report  the  functional  stress  intensity  factor  approach  is  summarized  for  both 
internally  cracked  and  externally  cracked,  partially  autof rettaged , 
pressurized  thick-walled  cylinders. 
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RESIDUAL  STRESSES  AND  THERMAL  SIMULATION 

The  residual  stress  distribution  in  an  autofrettaged  thick-walled 
cylinder  has  been  studied  by  a  large  number  of  investigators.  There  is 
considerable  disagreement  in  their  results  due  to  different  assumptions  which 
must  be  made  in  order  to  make  the  problem  mathematically  tractable.  Detailed 
discussions  of  the  results  and  the  associated  assumptions  are  given  in 
References  14  and  21.  Under  the  combination  of  assumptions  that  the  material 
is  incompressible,  elastic-perfectly  plastic,  and  obeys  the  von  Mises’  yield 
criterion,  and  that  the  cylinder  is  under  the  condition  of  plane  strain,  the 
following  closed  form  solution  for  residual  stresses  is  obtained  for  an 
elastically  unloaded  cylinder  after  partial  autofrettage : 


°r(r)  “  ^ 


(1) 


p  <  r  <  b  (2) 


oe(r)  -  < 


P  <  r  <  b  (4) 


1  <  r  <  p 


(3) 


where  bore  radius  is  taken  as  unit  length,  b  is  the  outer  radius  of  the  hollow 
cylinder,  p  is  the  radius  of  the  elastic-plastic  interface  during  pressuriza¬ 
tion,  aQ  is  the  uniaxial  yield  stress  in  tension  and  compression,  and 


This  residual  stress  distribution  is  used  in  this  report  as  a  basis  to 
develop  a  method  to  compute  stress  intensity  factors  for  cracks  in  such  a 
stress  field.  It  has  been  shown  in  Reference  13  that  the  thermal  stresses  in 
the  cylinder  subjected  to  a  thermal  load 


T(r) 


To  ~ 


(T0-Tp) 
log  p 


log  r 


l  <  r  <  p 


p  <  r  <  b 


(6) 


(7) 


are  equivalent  to  the  residual  stresses  In  Eqs.  (1)  through  (4)  If  the 
temperature  gradient  and  the  yield  stress  satisfy 

Ea(T0— Tp)  2 0O 

2(l-v)log  p  /3 

where  T0  and  Tp  are  the  temperatures  at  the  bore,  and  r  *  p  respectively,  E  is 
Young's  modulus,  and  a  is  the  coefficient  of  linear  thermal  expansion.  This 
thermal  simulation  provides  an  effective  method  to  compute  stress  intensity 
factors  due  to  Initial  stresses  given  by  Eqs.  (1)  through  (4)  using  NASTRAN 
or  APES.  The  initial  stress  is  replaced  by  a  temperature  input  of  Eq.  (6)  at 
all  nodes.  The  stress  intensity  factors  obtained  from  NASTRAN  or  APES 
corresponding  to  the  thermal  loads  are  equivalent  to  stress  intensity  factors 
due  to  autofrettage  residual  stresses. 22 


WEIGHT  FUNCTION  AND  FUNCTIONAL  STRESS  INTENSITY 

A  weight  function  is  a  universal  function  which  depends  only  on  geometry 
and  not  on  loadings. 20  If  the  mode  I  stress  intensity  factor  and 

displacement  field  associated  with  the  symmetric  load  system  1  are  known, 

the  weight  function  for  the  cracked  geometry  is 
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3u(1)(c) 


(8) 


2K( ^ (c)  8c 

where  H  -  E  for  plane  stress  and  H  ■  E/(l-v2)  for  plane  strain;  c  is  the  crack 

depth.  Once  h  is  determined,  the  mode  I  stress  intensity  factor  induced  by 

any  other  symmetric  load  system  x  and  f  is  given  by 

K  -  /  (T-h)dT  +  /  (f.h)dA  (8a) 

T  ~  ~  A 


where  t  is  the  stress  vector  acting  on  boundary  T  around  the  crack  tip  and  f 
is  the  body  force  in  region  A  defined  by  f.  This  equation  can  be  reduced  to 


H  c  3v(l) 

K  "  TiT  ^  Pc<x)  1 —  dx 

k(1)  O  3c 


(9) 


for  radially  cracked  rings  with  x  being  a  distance  measured  along  the  crack 
from  the  base  toward  the  tip.  The  relation  between  r  and  x  is 


1  4*  x 


for  interior  cracks 


r(x) 


(10) 


b  -  x  ,  for  exterior  cracks 

The  crack  pressure  pc(x)  can  be  found  from  the  hoop  stress  (at  the  site  of 
radial  cracks)  in  an  uncracked  ring  subjected  to  the  loading  of  interest. 

Even  though  the  numerical  values  of  and  v(l),  the  normal  component  of 

displacement,  are  known,  the  partial  derivative  3v(l)/3c  is  usually  unknown. 

A  technique  of  computing  3v/3c  was  devised  in  Reference  3  by  assuming  the 
crack  face  displacement  v  be  a  conic  section  given  by  Orange. 23  Another 
method  developed  in  Reference  15  made  no  assumptions  on  v  or  3v/3c,  but 
utilized  the  finite  element  method  to  compute  several  stress  intensity  factors 
each  associated  with  a  simple  loading  system.  For  a  new  load,  the  new  K  is 
expressed  in  terras  of  known  values  of  K. 
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The  hoop  stress  in  an  uncracked  cylinder  subjected  to  an  internal 


\ 

I 

I 


pressure  p^  is 


°6(r)  l  b2 

-  «  -r —  (1  +  -r) 

Pi  b2-l  r2 

Substituting  OQ  from  Eq.  (II)  as  pc  in  Eq.  (9)  we  have 


K(pi)  i 


—  Kc(l)  +  — 


Pi 


b2-l 


b2-I 


Kc(r"2) 


where 


H  c  3vU) 

M  1)  -  —  /  dx 

k( i)  °  3c 


H 


KcCr"2)  -  /  [r(x)] 


-2 


9vC  1 ) 


k(D  o  3c 

are  called  functional  stress  intensity  factors. 
Similarly  we  get 


dx 


K(p0) 


Po 


Kc(l)  + 

b2-l  b2-l 


(11) 


(12) 


(13) 


(14) 


(15) 


K(p-b)  1  0 

- -  —  U2  -  Pi(b)]Kc(l)  -  Pi(b)Kc(r-2)  +  2Kc(log  r)}  (16) 

o0  /3 


for  the  same  cylinder  subjected  to  uniform  tension  pQ  on  outer  cylindrical 


surface  and  fully  autof rettaged  residual  stress  respectively.  In  Eq.  (16)  the 


new  functional  stress  intensity  factor  is 


H  c  3v(D 

^(log  r)  -  ----  /  log( r(x) )  - - dx  (17) 

k( 1 )  o  3x 
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The  finite  element  results  of  K(pi)/pi,  K(p0)/p0  and  K(pab)/o0  enable  us 
to  compute  the  functional  stress  intensity  factors  1^(1),  K^r"2)  ,  and  K^log 
r) .  For  the  same  flawed  cylinder  with  a  different  degree  of  autof rettage ,  the 
stressses  can  be  computed  from  one  of  the  following  algebraic  equations.  For 
an  inner  crack  with  the  crack  tip  rc  in  the  range  1  <  rc  <  p,  the  equation  is 

K(p)  1  9 

-  =  —  {[2  -  Pi(p)]Kc(l)  -  Pi(p)Kc(r“2)  +  2Kc(log  r) }  (18) 

°o  /3 

For  an  outer  crack  with  rc  in  the  range  p  <  rc  <  b,  the  equation  is 

K(p)  1  ,  ,  . 

-  =  -  [p2-p1(p)j[b-2Kc(1)  +  KcCr"2)]  (19) 

°o  /3 


MODIFICATION  FORMULAS 

For  a  partially  autof rettaged  cylinder,  let  e  be  the  degree  of  auto- 
f rettage,  then  e  -  (p-l)/t  where  t  -  b  -1  is  the  wall  thickness  of  the 
cylinder.  When  a  crack  crosses  the  elastic-plastic  interface  r  *  p,  the  hoop 
stress  along  the  crack  face  must  be  represented  by  both  Eqs.  (3)  and  (4). 

Hence  Eqs.  (18)  and  (19)  are  not  valid  in  such  a  situation.  For  inner  cracks, 
we  may  use  Eq.  (18)  to  compute  an  approximate  value  which  is  based  on  the 
crack  face  loading  of  Eq.  (3).  The  error  introduced  by  Eq.  (18)  may  be 
corrected  by  adding  the  following  crack  face  loading 


Pc(r) 


0, 


1  <  r  <  p 


°o  r  p 2  1  1 

—  {2  log  -  +  1  +  -r  -  Pj(-r  +  -r)},  p  <  r  <  rc 
n  P  tr  r* 

(20) 


7 


where  rc  Is  the  radius  of  the  crack  tip.  Substituting  from  the  above  into  Eq. 
(9) ,  we  should  be  able  to  obtain  a  correction  stress  intensity  factor  K$  if 
3v/3c  is  known.  Assume  the  crack  tip  crosses  the  elastic-plastic  interface 
only  slightly,  the  Westergaard  near  field  solution  for  v  in  terms  of  crack-tip 
stress  intensity  factor  k(D  can  be  approximately  used  to  find  3v/3c.  Let  £ 
be  a  length  measured  from  the  crack  tip  and  be  defined  by 

5  -  -(x-c)  (21) 

then 


and 


2KU)  2  , 

v(0  - - (-)1/2  JT 

H  ir 


3v  kO)  2w.  1*1 

—  -  —  c-)1/2  (—  +  -  a) 

3c  H  TT  /£  C 


(22) 


(23) 


The  approximate  correction  factor  K$  obtained  by  termwise  integration  of  Eq. 
(9)  using  Eqs.  (20)  and  (23)  is  given  by 
Kfi  i 

—  -  --  fTR  (-1  +  2  log  p)(Ii  +  If)  +  pz(l2  +  12')  -  2(13  +  I3')>  (24) 

°o  /3 

Using  5t  ■  |rc-p|,  the  abbreviations  in  Eq.  (24)  are: 


II  **  2/ 6t  ,  lx'  “  —  (6t)^/2 

3c 

1  ^  (l+c)'1/2 

I2  =  -  [ - log  D(p)  ) 

1+c  p  2 


1 


^  (l+c)~1/2 


I2'  -  -  [ - + - log  D(p)  ] 

c  p  2 


(25) 


(26) 
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1+c  -  /St 


1 3  -  -2[(  2-log  p)/It  +  ( 1+c)  1/2  i0g - ] 

/1+c  +  /It 

(27) 

2 

I3*  -  —  [/It  { 6 1  log  p  -  2( 1+c)  -  26t/ 3}  -  ( 1+c) 3/2  log  D(p)) 

3c 

with 

D(p)  -  [2( 1+c)  -  p  -  2/ 6t( 1+c) )/p  (28) 

The  sura  of  K$/o0  from  Eq .  (24)  and  K(p)/o0  from  Eq.  (18)  usually  gives  a 
better  approximation  of  crack-tip  stress  intensity  factors  when  6  »  |rc-p|/t 
is  small.  Since  the  correction  formula  K^/o0  is  not  a  function  of  N,  the 
number  of  cracks,  it  vrorks  for  small  N  and  6.  But  when  N  is  large,  the  crack 
interaction  is  strong,  and  6  must  be  small. 

A  similar  formula  can  be  found  for  exterior  cracks  crossing  the  elastic- 
plastic  interface  from  elastic  into  plastic  region. 

2  , 

_  =  (__)1/2{(1_2  log  p)(Ji+Jl')  -  p2(J2+J2’)  +  2(J3+J3')}  (29) 

a  3tt 
o 

where 


2 

Jl  -  2/It  ,  Ji'  -  —  (6t)3/2 

3c 

(30) 

i  /£7 

J2  “  —  I - +  r^1/2  tan”V5t/rc] 

(31) 

rc  P 

i  /fiT 

J?'  “  “  l~  -  +  rc“l/2  tan-1/6t/rc] 

(32) 

c  p 

J3  *■  2[(-2+log  p)/It  +  2/r^  tan“V6t/rc] 

(33) 

2  2 

j3.  „  __  [(6t)3/2iog  p  +  -  (3rc-6t)/IF  -  2rc3/2  tan-Vet/rd  (34) 
3c  3 
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Adding  K$/0O  of  Eq.  (29)  to  K(p)/o0  of  Eq.  (19),  the  result  is  the  corrected 
crack-tip  stress  intensity  factors  for  exterior  cracks. 

NUMERICAL  RESULTS  AND  CONCLUSIONS 

Extensive  numerical  results  are  obtained  for  a  cylinder  of  b  *  2,  which 
is  a  commonly  used  value  in  cannon  design.  A  typical  finite  element  idealiza¬ 
tion  is  shown  in  Figure  1  for  inner  radial  cracks.  Slight  modifications  in 
element  meshes  are  needed  for  exterior  cracks.  In  Figure  1  the  elements 
surrounding  the  crack  tip  are  enriched  elements. If  collapsed  singular 
elements  are  to  be  used,  we  can  simply  replace  the  enriched  quadrilaterals  by 
collapsed  quadrilateral  elements  (triangular  elements)  with  proper  shifting  of 
side  nodes  to  new  locations.  Similar  accuracy  is  achieved  using  either 
enriched  or  collapsed  quadrilaterals.  Stress  intensity  factors  for  internal 
cracks  are  given  in  Table  1  of  Reference  15  for  three  different  types  of  load¬ 
ings.  Using  these  results  and  using  Eqs.  (12),  (15),  and  (16),  we  obtain 
functional  stress  intensity  factors  for  internal  cracks.  Figure  2  shows 
stress  intensity  factors  as  a  function  of  c/t  for  various  numbers  of  internal 
cracks  in  a  fully  autof rettaged  cylinder.  Figures  3,  4,  and  5  are  similar 
graphs  of  functional  stress  intensity  factors.  Corresponding  graphs  for 
external  cracks  are  shown  in  Figures  6  through  9. 

Readings  taken  from  Figures  3,  4,  and  5  are  enough  for  an  estimate  of 
stress  intensity  factor  for  internally  cracked  cylinders  with  any  assigned 
values  of  N  and  c/t  for  any  combination  of  p^  and  residual  stresses 
corresponding  to  a  given  e.  If  crack  tips  cross  the  elastic-plastic 
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interface,  then  correction  formula  (24)  should  be  used.  As  an  example,  if  the 
stress  intensity  factor  is  desired  for  N  -  2,  c/t  -  0,3  in  a  25  percent 
autof rettaged  cylinder,  we  first  take  readings:  K^O/Ziic  =*  1.41  from  Figure 
3,  KcCr"2)//^  -  1.05  from  Figure  4,  and  K^log  r)//rTc  “  0.22  from  Figure  5; 
then  KcCe-0. 25)/o0/iTc  *  -0.12  is  computed  from  Eq.  (18).  Since  rc  =*  1.3  is 
greater  than  p  *  1.25,  the  correction  stress  intensity  factor  K^/oq/ttc  * 

-0.023  is  obtained  from  Eq.  (24).  The  sum  of  Eqs.  (18)  and  (24)  gives  K/oq/ttc 
*  -0.143.  To  check  this  result,  a  finite  element  computation  of  this  case  is 
performed.  The  result  is  also  -0.143.  For  externally  cracked  cylinders,  Eq. 
(19)  involves  only  two  functional  stress  intensity  factors.  Therefore,  only 
Figures  7  and  8  are  needed.  For  example,  given  N  »  4,  readings  taken  from 
Figure  7  are  ¥^.(1) /Site  -  1.12  and  1.18  for  c/t  -  0.2  and  0.3  respectively. 
Readings  are  K^r*2)//*^  -  0.32  for  c/t  *  0.2  and  0.36  for  c/t  *  0.3  from 
Figure  8.  Stress  intensity  factors  for  z  ■  0.8  can  be  computed  from  Eq.  (19), 
giving  K(e*0.8)/oo/7Tc  -  0.492  and  0.537  for  c/t  *  0.2  and  0.3  respectively. 

For  c/t  *  0.3,  the  correction  formula  (29)  must  be  used.  The  result  is 
K^/oq/tTc  -  -0.05.  The  final  result  for  c/t  -  0.3  is  K(e«0.8)/oo  -  0.487  which 
is  close  to  the  result  of  0.486  obtained  directly  from  a  finite  element 
computation. 

For  a  combination  of  residual  stresses  and  internal  pressure,  the  stress 
intensity  factor  is  simply  an  algebraic  sum.  Stress  intensity  factors 
normalized  by  oQ/ Tc  are  shown  in  Figure  10  as  a  function  of  N  for  internal 
cracks  subjected  to  several  selected  values  of  p-t  and  e.  Figure  11  is  a 
similar  graph  for  external  cracks. 
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The  functional  stress  intensity  factors  are  used  to  obviate  the 
difficulty  in  finding  the  weight  function  itself.  The  approach  is  the  result 
of  a  series  of  methods  developed  for  multiply  cracked  cylinders.  The 
extension  of  the  method  to  residual  stress  distribution  other  than  that  given 
by  Eqs.  (1)  through  (4)  is  highly  possible. 

From  numerical  results,  the  stress  intensity  factor  is  largest  for  N  =  2 
for  various  combinations  of  residual  stresses  and  internal  pressures  for  both 
interior  and  exterior  cracks.  The  stress  intensity  factor  raonotonically 
decreases  as  the  number  of  cracks  increases  from  N  “  2. 
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Figure  1 (a) . 


(b) 


Figure  1(b). 


I 

CRACK  TIP 


A  typical  finite  element  idealization. 


CRACK  TIP 


Idealization  for  very  shallow  cracks 
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Figure  2.  Stress  intensity  factor  as  a  function  of  c/t  for  N  ID 
cracks  in  a  fully  autofrettaged  cylinder  of  b/a  =  2. 
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Figure  4. 
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Kc/pyric  vs.  c/t  for  N  internal  radial  cracks  with  crack 
face  loading  pc(x)  *  p(l+x)2> 
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Figure  5.  Kc/ p^ttc  vs.  c/t  for  N  internal  radial  cracks  with 
crack  face  loading  Pc(x)  ■  p  log(l+x). 
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Figure  6.  K/aQ^irc^  as  a  function  of  c/t  for  N  radial  cracks  at  outer 
surface  of  a  fully  autofrettaged  cylinder  of  b/a  =  2. 
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Figure  7.  Kc(p)/p/nx  as  a  function  of  c/t  for  N  external  radial  cracks 
with  constant  crack  face  loading. 
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figure  8.  Kc/p47  vs  c/t  for  N  external  radial  cracks  with  crack  face 
loading  ft  (x)  =  p(b-x)-2. 
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Figure  9. 


Kc/p/iiF  vs.  c/t  for  N  external  radial  cracks  with  crack  face 
loading  pc(x)  =p  log(b-x) . 
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Figure  11.  K/o0/iic  for  N  radial  cracks  at  outer  surface  of  a  cylinder  of 
b/a  =  2  subjected  to  combined  internal  pressure  pi  =  o0/f, 
where  f  *  1.5  except  otherwise  indicated,  and  residual 
stresses  corresponding  to  given  degrees  of  autofrettage  c. 
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